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A mechanical model for guided motion of mammalian cells 
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We introduce a generic, purely mechanical model for environment-sensitive motion of mammalian 
cells that is applicable to chemotaxis, haptotaxis, and durotaxis as modes of motility. It is able 
to theoretically explain all relevant experimental observations, in particular, the high efficiency of 
motion, the behavior on inhomogeneous substrates, and the fixation of the lagging pole during 
motion. Furthermore, our model predicts that efficiency of motion in following a gradient depends 
on cell geometry (with more elongated cells being more efficient). 

PACS numbers: 87.17.Jj, 87.17.Aa 


Motility and directed cell motion play an important 
role in many biological processes ranging from embry¬ 
onic development m to tissue invasion by pathogenic 
microorganisms [D13I1IS] and cancer progression [MS]- 

Often extracellular cues are used to regulate the deci¬ 
sion in which direction the cell will move |MII|. De¬ 
pending on these cues one distinguishes between: (1) 
chemotaxis where directed motion is guided by solvent 
chemical cues |9]; (2) haptotaxis where substrate-bound 
cues influence the cell-substrate adhesiveness (T^j; and 
(3) durotaxis where mechanical cues such as substrate 
rigidity influence the directed motion [Bj. 

Physics-based experimental and theoretical ap¬ 
proaches to study these phenomena have attracted 
considerable interest in the last decade. Many studies 
have been devoted to bacterial chemotaxis jasniiiMiT] 
and to Dictiostelium discoideum as model system for 
amoeboid migration [TTJ [T8l - t22] . The main challenge is 
to analyze and theoretically model the interplay between 
molecular processes and the emerging macroscopic mo¬ 
tion. For amoeboid motion this is further complicated 
as shape changes have to be taken into account pni422j . 

In contrast, the mesenchymal migration of mammalian 
cells has so far not been studied theoretically. Al¬ 
though many details have been characterized experimen¬ 
tally [IlESIEl] theoretical studies have focused only on 
the cell shape during migration |25|, on continuum de¬ 
scriptions for cell populations jZlIMj or on special short¬ 
term aspects of migration like migration speed |26j or 
effective adhesiveness [25] . 

In this paper we study the motility of mammalian 
epithelial-like cells. They are influenced by a variety of 
chemical and physical signals, in particular by different 
mechanical forces |27] and can show a very high efficiency 
in following all kinds of gradients [28]. Mammalian cells 
predominantly migrate by a crawling motion |29| . It con¬ 
sists of a cycle of five discrete steps carried out within 
about ten minutes [T1151151 1251 [2B1 [2BH55] : (1) polarization 
of the cell yielding a defined leading and a defined lagging 
pole; (2) formation of protrusions (predominantly) at 
the leading pole and attachment of the lagging pole, see 
Fig. Ha); (3) stabilization of these protrusions by adhe¬ 




FIG. 1: (a) Typical cell shape as observed in [^ for a cell po¬ 
larized in the direction of the arrow. The protrusions mainly 
grow into the forward direction, (b) In our model the po¬ 
larized cell is represented by an ellipse, the attachments by 
lines. The position of an attachment point is defined by the 
length of a protrusion and three angles: the rotation angle of 
the cell 9, the angle ip for the attachment position on the cell 
membrane, and the angle 7 between the central line and the 
arm. The distance between the center of the ellipse and the 
anchor point is 5. 

sion to the substrate or the extracellular matrix (ECM); 
(4) translocation of the cell-body by myosin-mediated 
contraction; (5) retraction of the rear by loosening the 
adhesions at the lagging pole. This crawling is propelled 
by the active lamellipodeum at the leading edge which 
pulls the passive cell body forward [M]. Cell motility 
depends on the stiffness of the substrate and the ECM. 
These effects are mediated indirectly via their impact 
on cell shape [55] and directly through so far unresolved 
mechanisms referred to as durotaxis [l[a[5M3H]. 

The mechanical interaction with the substrate or the 
ECM can be thought of as a bilinear sequential bind¬ 
ing [52], which affects the cell predominantly on the 
nanoscale through mechanosensing mechanisms [271 [55]. 
However, cells can also chemically manipulate the ECM, 
e.g., in case of cancer, where tumor cells stiffen the sur¬ 
rounding ECM [271 EZ] ) and build a rigid stroma around 
the tumor. This step in ECM stiffness then promotes 
cells from the outside moving inside, but prevents cells 
from the inside to migrate outside [55]. In general. 
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mammalian cells are not passive recipients of mechani¬ 
cal forces, but actively respond by pulling or pushing the 

ECM [SIETIISI]. 

We introduce here a simple, generic model for 
environment-sensitive motion of fibroblast-like cells. It 
is solely based on mechanics and is applicable to chemo- 
taxis, haptotaxis, and durotaxis as modes of motility. 
It provides the first theoretical explanation of the high 
efficiencies of mesenchymal-like motion independent of 
cellular morphology. Our model also covers the motil¬ 
ity dynamics on large time scales. In particular, we can 
capture the statistical properties at environmental dis¬ 
continuities, e.g., a step in substrate stiffness or a step 
in concentration of a chemo-attractant. These results 
indicate that regulation of taxis might be based on me¬ 
chanical forces. 

In our model we represent the polarized cell body by 
an ellipse with major radius Ri and minor radius Rg. The 
orientation with respect to the cc-axis is measured by the 
angle 6 that parameterizes rotation around the center C. 
To counterbalance the forces of the protrusions the cell 
is attached to the surface at an anchor point. Moving 
cells typically have a shape similar to the one shown in 
fig. [^a) with a single elongated tail that appears when 
the cell-body moves forward. This indicates that there is 
only a single anchor point. The motion occurs in such a 
way that the cell is effectively rotated around this anchor 
point, see Fig. m- In principle, cells could also fix the 
position of the tail with several anchor points. Then, 
the cell is effectively rotated around these fixed points. 
For simplicity, we restrict here the analysis to a single 
anchor point. We assume that the anchor point remains 
at a fixed position while the cell is not moving. Real cells 
might change their shape during motion that could lead 
to a change in the position of the anchor point. However, 
we do not take such effects into account. 

In the following we assume that a molecular mecha¬ 
nism initially polarizes the cell along the a;-axis (for ex¬ 
ample in the presence of a gradient as discussed below). 
For our purpose we do not explicitly model this process 
and assume that it leads to a normally distributed initial 
angle 0 (t = 0) with mean fig = 0. The standard devia¬ 
tion ag was estimated by fitting a normal distribution to 
the polarization model shown in [40) . 

The protrusions by which the cells pull themselves for¬ 
ward are represented by adhesive arms that grow out of 
the ellipse at a random angle ip with probability distri¬ 
butions given by Gaussian distributions centered around 
(f = 0 (for the leading pole) and (p = n (for the lagging 
pole) 

p{p) = ^ + Zr . (1) 

v 27 rCT+ v 27 rcr_ 

The weights of the leading pole and p_ = I —< p+ 
of the lagging pole reflect the initial polarization of the 


cell and shift the arm distribution towards the leading 
pole. 

Length and direction of the arms are random. For 
simplicity, we draw the angle 7 between growth direction 
and surface normal from a Gaussian distribution with 
mean /i-^ = 0 and standard deviation a^. Similarly, the 
arm length I is distributed normally with mean fii and 
standard deviation ai. 

Arm formation occurs at constant rate. Fvery arm 
applies a linear force on the cell that is proportional to 
the concentration of a chemoattractant or rigidity of the 
substrate. Thus, if at time t there are N arms that are 
attached to the cell body at position Xi = Xi{ipi) pointing 
in direction li (with I < i < N) the total force on the 
cell is given by 

N N 

F (rceii, 9,t) = '^F, = '^k (c) U (rceii, 9, Pi ), (2) 

i=l z=l 

where k (c) = kgc (rarm) depends on the chemoattractant 
concentration or the substrate rigidity c(r). This results 
in a translocation of the cell. Furthermore, the arms ex¬ 
ert a total torque M {rceii, 9, t) = ^ ^i)z^ 

leads to a rotation of the cell that can be interpreted as 
a gradual repolarization of the cell. 

For most of our simulations we used a linear gradient 
of fixed strength cq = 1/Ri- This results in a standard 
deviation of the initial angle of erg ~ 1 - 

To classify the planar cell motion (in the x-y plane), 
we numerically calculated the moments {x), {y),(x‘^), and 
( 2 /^), where (...) denotes the average over 500 indepen¬ 
dent runs, see Fig. In these simulations we imple¬ 
mented the above mentioned cycle of independent steps 
(l)-(5). Starting from the polarized shape [step (1)] we 
grow N new arms of lengths h at angles ipi in every itera¬ 
tion [steps (2) and (3)]. The new position and orientation 
of the cell is then obtained by solving F {rceii, 9, t) = 0 
and M {rceii, 9, t) — 0 independently [step (4)]. Then, all 
arms are removed [step (5)], the time is increased by 6t 
and the iteration starts over. In the absence of a gradient 
the cells perform an isotropic random walk. For iV = 10 
(for parameter values see [41]) the effective diffusion coef¬ 
ficients ^ Dy = 43.77i?;^/(5t were measured by fitting 
linear functions to (cc^) and (j/^). 

In this limit the model can also be solved analyti¬ 
cally. From the one-dimensional probability densities of 
ip, 7 and I the two-dimensional probability density of 
the force induced by an arm can be calculated [41]. We 
have compared these analytical results with those ob¬ 
tained by direct numerical integration of the model and 
found excellent agreement. It is interesting to note that 
this force probability density resembles the shapes of mi¬ 
grating lamellipodial domains of keratocytes [25]. If one 
assumes that these shape deformations reflect the forces 
then the forces acting on the cell in our simple proba¬ 
bilistic model are remarkably similar to those exerted by 
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FIG. 3: Efficiency of motion (as quantified by the chemotactic 
factor CF) as function of the number of arms N for different 
simulation time steps n. 


FIG. 2: Glassification of ceil motion as obtained from the nu¬ 
merically calculated moments {x), and (y^). In the 

absence of gradients the cell performs a random walk charac¬ 
terized by vanishing {x) and (y) and standard deviations that 
increase with time. In the presence of a gradient this random 
walk becomes biased. In this case the motion in y-direction is 
suppressed compared with the case of no gradient. Data were 
obtained by averaging over 500 independent runs. Length and 
time are given in units of Ri and 5t, respectively. 


the actin cytoskeleton of keratocytes. 

Next, we consider substrates with gradients. For a lin¬ 
ear gradient cq = l/Ri parallel to the x-axis, the cells 
perform a biased random walk in the gradient direction 
while the motion in the perpendicular direction is sup¬ 
pressed ^ 135-fold [Dy = 0.32R‘f/6t in presence of a 
gradient compared to Dy = i3.77Rf/St in absence of a 
gradient), see Fig. To quantify the efficiency of the 
motion in following the applied gradient we measured 
the chemotactic factor 


balance between increase of membrane area and gain of 
efficiency within this range of N. This number is also 
comparable to the number of arms found experimentally 

El- 

Next, we analyzed the influence of the cell geometry on 
the efficiency of motion. As Fig. shows, CF depends 
on the geometry of the cell characterized by the ratio 
Ri/Rs- Thus, more elongated cells (with Ri > Rs) have 
a higher CF. This is somewhat surprising as these cells 
have a broader force distribution than less elongated cells. 
However, as we show in m CF depends predominantly 
on the ability of the cell to align with the prescribed gra¬ 
dient. Thus, for more elongated cells this higher ability 
compensates for the broader force distribution. The CF 
as determined by the ability to align with the gradient 
(characterized by a rotational rate a) is given by [IT] 


CF = exp(x^) erf ^TT/{V2ag) — + erf(x) 


(4) 


CF = 


Agrad \ 

Rtot j 


(3) 


Here, Ltot is the total path length and Lgrad the length 
of the projection in the direction of the gradient. 

To investigate the robustness of our model to vary¬ 
ing gradients cq and its behavior for small gradients we 
looked at the dependence of CF on cq, see Fig. S4 in mi- 


For small gradients we see a strong increase in efficiency 
with the gradient strength, but the efficiency saturates 
fast to its maximum value of around CF = 95%. 

An increase in number of arms results in a speedup of 
motion and an increase in efficiency, see Fig.[^ However, 
there is saturation in efficiency and speed for large num¬ 
bers of arms. If we take into account that 5 to 10 arms 
with an average length of 7/im would roughly cover be¬ 
tween 10 and 20% of the surface of the cell (each arm has 
about 22/im^ surface area while the whole cell body has 
about lOOO^m membrane area [321133]) we reach a good 


Here, x = 2~^/'^aga/x and the opening angle of the 
force distribution y = arctan (i?//i?s tan (t+) parameter¬ 
izes the dependence on cell geometry. 

Lo et. al. |5] have shown that non-moving cells grow 
longer protrusions on stiff substrates than they do on 
soft substrates. This implies a regulative effect of the 
substrate rigidity. To account for this effect in our model 
we assign each arm i the length 


li -j- Sli (c) — li 


SL 


1 _l- ^ ■ 


(5) 


where the regulated elongation Sk (c) depends on the 
stiffness (or concentration or adhesiveness) c = c (rarm) 
at the position of the arm rami- 

If we increase the average arm length, we see an in¬ 
crease in efficiency that saturates for longer arms, see 


Fig. S5 in [41j . This rise in efficiency comes with an in¬ 
crease in speed. If we use the concentration-dependent 
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FIG. 4: Chemotactic factor CF is shown as function of cell 
geometry characterized by Ri/Rs- More elongated cells follow 
with higher efficiency the gradient (of strength cq = 1/Ri). 
The data shown are for different simulation time steps n. For 
smaller n (i.e. earlier simulation times), CF still depends on 
the initial alignment of the cells. With increasing n the cells 
align with the gradient increasing in this way CF. However, 
as can be seen only cells with x > 0.06 (i.e. cells which are 
sufficiently elongated) are able to align with the gradient. 


FIG. 5: Dependence of CF on the position of the anchor 
point to the substrate. The distance of this point from the 
center of the ellipse is denoted by 5. As this point moves from 
the center of the ellipse towards the leading pole, CF drops 
sharply to zero. This indicates that in order to be able to 
follow a prescribed gradient the cells need to have a defined 
leading pole and an anchor point sufficiently far away from 
this pole. 


regulation of arm lengths we see a 1 % increase in maxi¬ 
mum efficiency compared to the unregulated system, but 
the same efficiency is reached with an average arm length 
up to 40% shorter compared to the unregulated system. 

The above results are robust with respect to variations 
in the standard deviations of the distributions for ip and 
7 . CF remains nearly constant for (T+ and in a range 
from ~ O.OItt to 0.27r. For even broader distributions we 
see a decrease in efficiency as a result of insufficient polar¬ 
ization of the cell. Experimentally, it has been observed 
that the protrusions grow almost in the normal direction 
out of the surface of the cells, and that the protrusions 
are located around a narrow region at the leading pole 
0 ISM]- We find a similar behavior (and the asso¬ 
ciated high efficiencies in motion) only for distributions 
with small standard deviations indicating that ip and 7 
are tightly regulated. 

Furthermore, cell motion shows an interesting depen¬ 
dence on the position of the anchor point of the cell to 
the substrate quantified by the parameter 5. If this point 
is shifted towards the leading pole, i.e. closer to the pro¬ 
trusions at the front, the torque exerted by these arms is 
reduced due to the shorter lever arm and the efficiency 
drops sharply to zero indicating that the cell is not able 
to follow the gradient at all (see Fig. [^. This indicates 
that the efficiency in following a gradient is dominated 
by protrusions close to the leading pole. 

Cells often encounter inhomogeneous substrates. As 
a general scenario, we analyze the movement towards a 
step in substrate rigidity where crawling cells show an 
interesting behavior. At these steps (that could repre¬ 
sent the transition from a rigid stroma of a tumor to 


the softer surroundings M) cells tend to move from the 
softer substrate to the stiffer substrate [39]. Cells mov¬ 
ing in this direction are only weakly influenced by the 
step. They keep moving but their trajectory bends to¬ 
wards the direction perpendicular to the step, see Fig. [^ 
We observed that the relation between the angle of the 
cell before and after the step obeys a refraction law simi¬ 
lar to that of light allowing us to characterize the motion 
by refraction indices, see Fig. [^a). The ratio of refrac¬ 
tion indices u-soft/’T'Stiff = sin^out/ sin 0 i„ decreases with 
increasing step size s = cijcQ, where cq and ci are the 
stiffness of the softer and stiffer region, respectively. 

On the other hand, a step from a stiff substrate to a 
softer substrate represents a barrier for cells. The pass¬ 
ing probability depends on the step height. From the dis¬ 
tribution of the minimal ^-positions encountered by the 
cells during 500 iterations, we can calculate the proba¬ 
bility of a cell moving across the step, which we define as 
transmission coefficient Tc = p (x) = p[x < a:step), 
where p (x) is the probability to find a cell that trav¬ 
eled to position x, and a:step is the position of the step, 
see Fig. [^b). This coefficient decreases as the step size 
increases showing that the barrier effect becomes much 
stronger for larger steps, see Fig. S7 in mi- 


Finally, we wanted to compare our results with other 
types of cellular motion. However, there are neither 
theoretical nor experimental data available for the effi¬ 
ciency of mammalian cell motion as function of gradient 
strength. We therefore decided to compare our model 
with data for swim-tumble chemotaxis, the standard 
model for cellular motion. More specifically we tested our 
results against a two-dimensional swim-tumble model for 
bacterial chemotaxis with a parameter set optimized for 
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^min 


FIG. 6: Transition at a sharp step, (a) If cells move from the 
soft region to the stiff region one hnds a refractive behavior 
depending on the step size s. (b) Cells moving from the stiff 
to the soft region experience a barrier that they overcome 
with a step-size-dependent transmission coefficient. For large 
step sizes no cells cross the barrier while for small step sizes 
almost all cells are able to cross it. 
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Co 

FIG. 7: Variation of the efficiency with the gradient strength 
cq. The cells reach their maximum efficiency at rather small 
gradients compared to a swim-tumble model of bacterial 
chemotaxis (based on [9] [10]). 


efficiency [adoiiiT]. As Fig.j^shows our model yields sig¬ 
nificantly higher efficiencies than the swim-tumble model. 

To summarize, in this study we have introduced a 
generic mechanical model for environment-sensitive mo¬ 
tion of mammalian cells. Motion occurs by polarized 
growth of protrusions which push and rotate the cell. 
The description of molecular interactions occurs on a 


coarse-grained level effectively entering into the proba¬ 
bility distribution for arm growth, see Eq. ([^. It is not 
the goal of this simplified model to achieve a detailed 
(molecular) description of cellular motility. Rather, we 
introduce it to analyze the influence of mechanical forces 
on the regulation of fibroblast motion. 

The agreement of our findings with the experimental 
observations indicates that mechanical forces indeed play 
a significant role in this process. More specifically, the 
observed high efficiency in following a gradient is a ro¬ 
bust feature of our model (see figs. and [^. We find for 
a large parameter range chemotactic factors close to 1 
as observed experimentally in [28] . Furthermore, the re¬ 
sults on the motion in inhomogeneous environments are 
in good agreement with the experimental observations at 
steps in substrate rigidity [39] . 

Our analysis identifies two geometrical factors that 
have a significant impact on efficiency of cell motion: the 
position of the anchor point and the geometry of cell 
quantified by the ratio of major to minor radius. For 
both quantities we make specific theoretical predictions 
(figS-H and[^ that are experimentally testable. 

Furthermore, the moments (x), (x^), (y) and (y^) can 
be easily measured for individual cells of different geome¬ 
try for different gradients. By comparing these data with 


in |41| ) information can be obtained about the concen¬ 
tration dependent regulation of arm lengths. To check 
our results concerning steps in concentration, adhesive¬ 
ness and stiffness, one could use the methods presented 
in to produce flat substrates of different stiffness and 
measure the polarization axes of the cells before and af¬ 
ter the interface as well as the transmission coefficients 
with time-lapse microscopy. 

There are many possible extensions of our model. In 
future work we will take into account the mechanical ef¬ 
fects that the cells have on the substrate. If cells attach 
protrusions to the substrate and contract, they locally 
stiffen the substrate. This local stiffening of the sub¬ 
strate might lead to an effective attraction between two 
cells in proximity, in this way promoting aggregation. 


our theoretical predictions (Fig. m and Figs. S5 and S6 
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Supplemental Material: A mechanical model for guided motion of mammalian cells 


ANALYTICAL RESULTS 


Force distribution 


The cell is represented by an ellipse with semiaxis Ri > Rs (see Fig. SI I. An adhesive arm grows out of the ellipse 
at angle (p. We now choose a coordinate system (a;®®', with the fixed counterpart of the protrusions as origin and 
the semimajor axis as abscissa. The vector to the origin of an arbitrary arm is given by 


Ri (cos p — S) 
Rs sin If 


The angle between the vector to the attachment point and the abscissa is 

Rs sin(vj) 


ip' = arctan 


Ri {cosp — (5) 


(S.l) 


(S.2) 



FIG. SI: Ellipse in the coordinate system with fixed counterpart of the protrusions as origin and the semimajor axis as abscissa. 
Vectors to the attachment point from the origin and from the center of the ellipse are shown in green. The coordinate systems 
in the attachment point are shown in red: solid lines for the system with one axis perpendicular to the surface, dashed lines 
for the system rotated parallel to the system of the semiaxis. 


Length I and direction 7 of the arms are Gaussian distributed 

1 


p{j) = 




■ exp 



( 

2 



p{i) = 


V^cri 


exp 


I fl-l 

'2 


(S.3) 


Here, the angle 7 is measured relative to the normal of the ellipse in the attachment point. The probability to find 
an arm with length I and angle 7 is 


=p{l)p{l) = 




TT (X-T, 


• exp 



/ 

2 

_ 




exp 


1 / 1 -/ 

'2 


The force caused by an arm with length I and angle 7 is given by 


= kl 


cos 7 
sin 7 


P-L I ’ 

V 


(S.4) 


(S.5) 


in the coordinate system {x_\_,y±) in the attachment point (see Fig. SI). In the absence of a gradient we can assume 
k=l. 


























FIG. S2: Probability density for the force ( Fx, Fy ). (a) Numerical result; (b) analytical result (S.91 (for same set of parameters); 
(c) analytical result in the approximation (|S.10 1. 


The probability to have a force {F^,F^) is 
1 1 




2na,a, ^+ F^^^ 


exp 




2 n 




exp 


1 / arctan(F^;L,F^) 


where we denote by arctan(M,r;) the arc tangent of ^ taking into account the quadrant of the point {u,v). 
To get the force in the coordinate system of the semiaxes, we need to rotate by the angle 

Ip = arctan (i?/ sin ip, Rg cos (p). 

Thus, 

p^,^{F,x,Fy) = p-^{F^,Fy) = p-^{Fx, cos'ip + Fy s\w'ip,-Fx s\n'ip + Fy cosip) 


,Jf^ + F 2 


exp 


1 I \/F S + FS- l 
m 


X exp 


1 / arctan(—Fj; Ri sin p + Fy Rg cosp,Fj; Rg cos p + Fy Ri sin(/3) 

2 V a-, 


This probability density is valid for any (fixed) p. The complete density is given by 

TT 

P-(F„ Fy) = J dp P^;iF,, Fy) pip), 


(S. 6 ) 


(S.7) 


(S. 8 ) 


(S.9) 


which has to be calculated numerically. The result shown in Fig. |S2| is in good agreement with the numerical data. 
As Fig. S2 shows the distributions of 7 and p are rather sharp (i.e. a^,- < ^ 1). Therefore, we can approximate 


the corresponding Gaussian functions by delta functions leading to 
f’^iF^.Fy) 


RiRg^F^ + F^ 


27ra,a^+ R^F^ + RlF^ 


■ exp 


1 Nf^ + f^-1 


<Ji 


exp 


1 ^arctan (i?s Fy,Ri F^,) 

'2 


<Jip+ 


+ 


P- 


1 


Jf^ + F 2 


exp 



exp 


1 /arctan (—Fy,—Fa;) 

2 [ 0-7 '^ 


(S.IO) 
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From this approximation, we conclude that the width of the distribution is given by ct/ in radial direction and by 
and cr-y in angular direction. We define x the angle, for which arctan(i?s Fy,Ri F^) < cri^+, i.e. 

X = arctan (—- tancr,p_|_ ) . (S-H) 

KRs ) 

Thus, X quantifies the ability of the cell to move straight ahead and/or turn towards a given gradient. It depends on 
cr,p+ and on the geometric factor 


Model for motion along gradients 

To investigate the influence of the opening angle on the ability to follow a gradient we analyze a somewhat simplified 
scenario: we neglect all effects from the back of the cell and do not take into account the torque. 

In a first step we describe the motion as a random walk. In doing so we keep track of the position and direction 
of cell. In each step the cell follows its direction for a constant distance and then changes direction by an angle ±77 
relative to its current direction. The probabilities for selecting +ry and —77 are denoted by p and 1 — p, respectively. 
The probability that after N steps the cell is rotated n times by 77 and N — n times by —77 (resulting in a net rotation 
of {2n — N)r]) is given by 


P7v(n)= (^) p" (S.12) 

For large numbers {N ^ 1, tt 1, fV — tt 1) this becomes a continuous Gaussian distribution for the net rotation 
X = (2n — N)r] 


P{x) 


,— exp 

\f^a 


1 {x — x)^ 

2 


(S.13) 


with mean x = N{2p — 1)77 and standard deviation a = 2^J A^p(l — p)r]. 

The efficiency of the motion can be quantified by the expectation value of cos x which represents the fraction of the 
distance covered in the direction of a gradient in a;-direction 


(cos x) 


= dx cosx 




■ exp 




(S.14) 


In the last equation one can identify the standard deviation cr with the opening angle x of the force distribution. Thus, 
this simple model predicts that with increasing x the efficiency should decrease. As can be seen from Fig. |S3| this 
clearly contradicts our numerical findings. This indicates that the ability to follow the gradient is not determining the 
efficiency of motion. As we show now, it in fact depends crucially on the ability of the cell to align with the gradient. 

To do so, we assume that initially the cell has an angle 9q between its major axis and the applied linear gradient. 
With a probability p the cell now rotates by an angle 'k towards the gradient. Thus, after n = ^ the cell is perfectly 
aligned with the gradient. 

Here, the efficiency can be quantified by the averaged probability of the cells to align with the gradient 


TT 



— TT 


(S.15) 


Taking into account that the starting angles 9q are Gaussian distributed in our simulations, the efficiency becomes 


d9c 




1 9. 


IM 


= exp 


2 af 

r gg (Inp) ^l 
2 ^-2 


y/2'Kag 


erf 


1 0 ; 


exp ( ^ ) pA' = ——— / f70o exp ( ^ + “t 


2 al • ^ 


V^crg 


o-g Inp 
4- 


erf 


<Jg Inp 

V 2 4- 


(S.16) 
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X 


FIG. S3: Efficiency of motion as function of x a-nd number of time steps n in the simulation. At early simulation stages 
the efficiency is dominated by the alignment with the gradient. At later times all cells are aligned with the gradient and the 
efficiency approaches 1. For vertically elongated ellipses {Rs > 2Ri) the torque can turn the cell against the gradient resulting 
in a negative efficiency. If this event occurred the simulation was stopped. The lines show the results of htting Eq. (S.16I to 
the numerical data. 


Thus, the probability depends on the rotation angle T* (that we identify with the opening angle x of tfi® force 
distribution) and on the probability p (which depends on the strength of the gradient). Because the exact relation 
between p and the strength of the gradient is not known we take Inp as a fitting parameter. As one can see from 
Fig. |S3| the fitted curves match the numerical data quite well for sufficiently large number of simulation steps. At 
early simulation stages the numerical data depends still on the initial orientation of the cell. 

For vertically elongated cell {Rs > 2Ri) the gradient eventually turns the cell into the “wrong” direction. This 
tilt in the wrong direction increases with time since arms growing in this direction are favored by the gradient. This 
process leads to an alignment of the ellipse against the gradient. We stopped our simulations when such an event 
occurred. 


MODEL FOR BACTERIAL CHEMOTAXIS 

To model the bacterial chemotaxis, we use the well-studied swim-tumble model. The cell swims into a random 
direction for a normally distributed length 1. It measures the change in concentration Ac of an attractant along this 
path. On the basis of Ac the decision is made whether the cell maintains the current direction or tumbles. The former 
occurs with probability 


Po 


Ac(l -po) 


1 - 


p{Ac) = < 


1 


kpo 


1 + k{l-po) 

Ac -I- Po (k — Ac) 


Ac < 0 


Ac > 0 


(S.17) 


which is larger than a basal probability of p (Ac = 0) = po for Ac > 0 and smaller than po for Ac < 0. The parameter 
k quantifies the stiffness of the response. 

To compare this model with our model, we also applied here a concentration gradient with a linear slope in x- 
direction. 
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SUPPLEMENTARY FIGURES 



FIG. S4: Gains in chemotactic factor CF, i.e. difference between CF at the end and at the beginning of the simulation, as 
function of gradient strength co. Results are shown for different standard deviations ag of the initial polarization angle Oq. For 
smaller ug the cells are initially aligned more in the direction of the gradient, so less gain is expected. 



eff ' i*! eft 


FIG. S5: Effect of the concentration dependent regulation of arm lengths on efficiency of motion, (a) If we increase the 
average arm length fii we see an increase in efficiency that saturates for longer arms. If we use the concentration dependent 
regulation of arm lengths see Eq. we get a 1% increase in maximum efficiency compared to the unregulated system 

from CF ~ 0.95 to CF ~ 0.96, but the same efficiency is reached with an average arm length up to 40% shorter compared to 
the unregulated system, (b) Speed of motion, characterized by the covered path length Ax per time unit At increases with 
increasing arm length. Additionally, for identical average arm length motion is faster with regulation. The difference in speed 
increases with increasing average arm length. 
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t/io'* 


FIG. S6: Classification of the migration dynamics. The cell performs a random walk which is biased in the presence of a 
gradient, see Fig. If the arms are regnlated the motion in absence of a gradient stays the same, whereas the motion in 
presence of a gradient becomes faster at constant efficiency (additional red dash-dotted line). 




FIG. S7: Histogram and probability of intrusion at a surface step (i.e. interface between stiff and soft substrate region). Cells 
move from left (grey, stiff substrate) to right (white, soft substrate). The histogram shows the frequency of the maximal (most 
right) X value reached. We only show the region close to the step to focus on the cells that did not cross it. Cells that have 
passed the step (most of the cells in simulations for small steps, i.e. blue and green) would be found at Xmax > 2, so they 
are not shown in the histogram. For larger step sizes not all cells pass the step. The maximal distance traveled in positive 
a;-direction peaks closer to the interface with decreasing step size. If the step size becomes too small almost all cells pass the 
step as indicated by the intrusion probability that is defined as the probability that a cell travels further than Xmax. 
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SUPPLEMENTARY TABLES 


Parameter Fig. |2 & S 6 Fig. 3 Fig. 


& 


S3 


Fig. 5 Fig. 6 'a) Fig. 6 'b) Fig. 7 Fig. S2 Fig. S4 Fig. S 6 Fig. 


S7 


Iter. 

81000 

- 

- 

500 

500 

500 

500 

1 

50C) 

500 

500 

Runs 

500 

500 

1000 

500 

500 

500 

500 

5000000 

500 

500 

500 

a 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

b 

0.5 

0.5 

- 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

5 

-1 

-1 

-0.9 

- 

-1 

-1 

-1 

-0.9 

-1 

-1 

-1 

CTO 

1.047 

1.047 

1.047 

- 

- 

0 

1.047 

- 

- 

1.047 

0 

k 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Co 

- 

1 

1 

1 

1 

1 

1 

- 

- 

1 

1 

N 

10 

- 

10 

10 

10 

10 

10 

1 

10 

10 

10 

P+ 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

Pi 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

- 

0.5 

(71 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 


0.5 

0.1 

0.1 

0.5 

0.5 

0.5 

0.5 

0.3 

0.5 

0.5 

0.5 

G- 

0.05 

0.01 

0.01 

0.05 

0.05 

0.05 

0.05 

0.03 

0.05 

0.05 

0.05 


0.7 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 


TABLE SI: Parameter values for the results shown in Figs. 2-7, S2-S7 


Step size s 

^Soft 2^- 

2 

0.9374 1.0000 

5 

0.8843 0.9480 

10 

0.8622 0.3220 

100 

0.8330 0.0020 


TABLE S2: Values for the ratio of refraction indices ttSsstt and transmission coefficient Tn. The refractive effect becomes 
stronger with increasing step size whereas the transmission decreases. 













































